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AIRBORNE WIND PROFILING ALGORITHM 
FOR DOPPLER WIND LIDAR 

CROSS-REFERENCE TO RELATED PATENT 

APPLICATIONS 5 

This patent application claims the benefit of and priority to 
U.S. Provisional Patent Application Ser. No. 61/751,291, 
filed on Jan. 11, 2013, U.S. Provisional Patent Application 
Ser. No. 61/809,664, filed on Apr. 8, 2013, U.S. Provisional 
Patent Application Ser. No. 61/820,478, filed on May 7, 2013, 
and U.S. Provisional Patent Application Ser. No. 61/820,957, 
filed on May 8, 2013. This application is also related to 
co-pending U.S. patent application Ser. No. 14/079,965, 
titled “Airborne Doppler Wind Lidar Post Data Processing 
Software DAPS-LV,” filed on Nov. 14, 2013. The foregoing 
applications are hereby incorporated by reference in their 
entirety. 
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STATEMENT REGARDING FEDERALLY 
SPONSORED RESEARCH OR DEVELOPMENT 

The invention described herein was made by employees of 
the United States Government and may be manufactured and 25 
used by or for the Government of the United States of 
America for governmental purposes without the payment of 
any royalties thereon or therefore. 

BACKGROUND OF THE INVENTION 30 

Light Detection and Ranging (“LIDAR”) is a sensing tech- 
nology using light pulses (e.g., light from a laser) to deter- 
mine a distance to objects. Doppler Wind LIDAR systems use 
light pulses to generate measurements of wind fields in the 35 
atmosphere based on the characteristics of the reflected light 
pulses returned from aerosols suspended in the atmosphere. 

For example, Doppler Wind LIDAR systems may calculate 
the distance to aerosols suspended in the atmosphere based on 
the time difference between transmitting a light pulse or 40 
signal and receiving the reflection of that light pulse or signal 
and/ or may calculate the velocity of the aerosols suspended in 
the atmosphere based on the frequency shift between a trans- 
mitted light pulse or signal and a received reflected light pulse 
or signal. The measurements generated by Doppler wind 45 
LIDAR systems may be used to generate wind profiles and 
other models providing information about the atmosphere. 

Airborne Doppler wind LIDAR systems (e.g., Doppler 
wind LIDAR systems operating from aircraft, such as planes, 
moving through the atmosphere) present challenges in 50 
extracting wind parameters due to many factors such as: data 
including the movement of instruments (e.g., the instruments 
to transmit and receive light pulses move with the aircraft and 
the data is impacted by the instrument movement); wind 
parameters (e.g., wind speeds) being smaller in magnimde 55 
than motion parameters of an aircraft (e.g., speed, roll, pitch, 
and yaw distortion); and instruments not being perfectly 
mounted in the aircraft. In order to accurately extract wind 
parameters from airborne Doppler wind LIDAR system data, 
these factors must be accounted for in an airborne wind pro- 60 
filing algorithm. 

BRIEF SUMMARY OF THE INVENTION 

The systems, methods, and devices of the present invention 65 
enable airborne Doppler Wind LIDAR system measurements 
and Inertial Navigation System/Global Positioning System 
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(“INS/GPS”) measurements to be combined to estimate wind 
parameters and compensate for instrument misalignment. 

The systems, methods, and devices of the various embodi- 
ments provide an Airborne Wind Profiling Algorithm for 
Doppler Wind LIDAR (“APOLO”), which is a breakthrough 
in airborne wind profiling for Doppler Wind LIDAR systems. 
In an embodiment, APOLO may have two components, a 
wind parameter estimation component and an instrument 
misalignment compensation component. In another embodi- 
ment, APOLO may estimate five instrument installation bias 
angles such as INS/GPS roll, pitch, and yaw biases, as well as 
azimuth and elevation biases. In yet another embodiment, 
APOLO may accurately estimate horizontal wind speed and 
direction at different altitudes for airborne Doppler Wind 
LIDAR systems. In an alternative embodiment, APOLO may 
extract wind parameters from other airborne wind LIDAR 
systems as well. 

In yet another embodiment, APOLO may function by (I) 
estimating the instrument installation bias angles such as 
INS/GPS roll, pitch, and yaw, and elevation and azimuth 
biases; (2) given an airborne Doppler Wind LIDAR data and 
an INS/GPS data file, distributing the data in range gate bins; 
(3) estimating the power spectra of the monitoring pulse data 
in the beginning of each wind LIDAR data; (4) performing 
DVT (“Data Validity Test”) by investigating the peak fre- 
quency of the monitoring pulse data; (5) estimating the power 
spectra and the Doppler shift of each range gate bin data using 
their periodogram; (6) transforming the wind vector from the 
Earth coordinate (East -North-Up) into the aircraft coordinate 
using the INS/GPS data and the estimated instrument instal- 
lation bias angles; (7) estimating wind parameters in the 
aircraft coordinate; and (8) transforming the wind parameters 
back to the Earth coordinate. 

In a further embodiment, the wind speed and wind direc- 
tion may be computed based on two orthogonal line of sight 
LIDAR remms. In another embodiment, the wind speed and 
wind direction may be computed based on five line of sight 
LIDAR returns. 

These and other features, advantages, and objects of the 
present invention will be further understood and appreciated 
by those skilled in the art by reference to the following speci- 
fication, claims, and appended drawings. 

BRIEF DESCRIPTION OF THE SEVERAL 
VIEWS OF THE DRAWINGS 

The accompanying drawings, which are incorporated 
herein and constitute part of this specification, illustrate 
exemplary embodiments of the invention, and together with 
the general description given above and the detailed descrip- 
tion given below, serve to explain the features of the inven- 
tion. 

FIG. 1 is system block diagram of an example airborne 
Doppler Wind LIDAR system suitable for use with the vari- 
ous embodiments; 

FIG. 2 is a component block diagram of an example data 
acquisition processing and display system according to an 
embodiment; 

FIG. 3 is a data structure diagram of example modules of 
data acquisition and processing software according to an 
embodiment; 

FIG. 4 is a process and data flow diagram illustrating opera- 
tions performed by an airborne Doppler Wind LIDAR system 
during airborne wind profiling according to an embodiment; 

FIGS. 5A-5D are a process flow diagram illustrating an 
embodiment method for computing wind speed and wind 
direction from Airborne Wind LIDAR data; 
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FIGS. 6A-6D illustrate various relationships between vari- 
ous aircraft and wind components, such as speed, heading, 
pitch, roll, yaw, and direction and aircraft and Earth coordi- 
nates; 

FIG. 7 illustrates a comparison between graphs of APOLO 
wind speed data and dropsonde wind speed data before and 
after applying offset compensation; 

FIG. 8 illustrates a comparison between graphs of APOLO 
wind direction data and dropsonde wind direction data before 
and after applying offset compensation; 

FIG. 9 illustrates a comparison between graphs of 5 LOS 
based wind speed data and dropsonde wind speed data using 
different offset estimates; 

FIG. 10 illustrates a comparison between graphs of 5 LOS 
based wind direction data and dropsonde wind direction data 
using different offset estimates; 

FIG. 11 illustrates a comparison between graphs of 
APOLO and a 5 LOS based wind speed data and dropsonde 
wind speed data with offset compensation applied; 

FIG. 12 illustrates a comparison between graphs of 
APOLO and a 5 LOS based wind direction data and drop- 
sonde wind direction data with offset compensation applied; 

FIG. 13 is a graph of aerosol layer distribution over time 
based on APOLO data witli offset compensation applied; 

FIG. 14 is a graph of wind speed over time based on 
APOLO data with offset compensation applied; and 

FIG. 15 is a graph of wind direction over time based on 
APOLO data with offset compensation applied. 

FIG. 16 is a component block diagram of a computing 
device suitable for use with the various embodiments. 

DETAILED DESCRIPTION OF THE INVENTION 

For purposes of description herein, it is to be understood 
that the specific devices and processes illustrated in the 
attached drawings, and described in the following specifica- 
tion, are simply exemplary embodiments of the inventive 
concepts defined in the appended claims. Hence, specific 
dimensions and other physical characteristics relating to the 
embodiments disclosed herein are not to be considered as 
limiting, unless the claims expressly state otherwise. 

The word “exemplary” is used herein to mean “serving as 
an example, instance, or illustration.” Any implementation 
described herein as “exemplary” is not necessarily to be con- 
strued as preferred or advantageous over other implementa- 
tions. 

The various embodiments will be described in detail with 
reference to the accompanying drawings. Wherever possible, 
the same reference numbers will be used throughout the 
drawings to refer to the same or like parts . References made to 
particular examples and implementations are for illustrative 
purposes, and are not intended to limit the scope of the inven- 
tion or the claims. 

The systems, methods, and devices of the present invention 
enable airborne Doppler Wind LIDAR system measurements 
and INS/GPS measurements to be combined to estimate wind 
parameters and compensate for instrument misalignment. 

The systems, methods, and devices of the various embodi- 
ments provide an Airborne Wind Profiling Algorithm for 
Doppler Wind LIDAR (“APOLO”), which is a breakthrough 
in airborne wind profiling for Doppler Wind LIDAR systems. 
In an embodiment, APOLO may have two components, a 
wind parameter estimation component and an instrument 
misalignment compensation component. In an embodiment, 
APOLO may estimate five instrument installation bias angles 
such as INS/GPS roll, pitch, and yaw biases, as well as azi- 
muth and elevation biases. In an embodiment, APOLO may 


4 

accurately estimate horizontal wind speed and direction at 
different altitudes for airborne Doppler Wind LIDAR sys- 
tems. In an embodiment, APOLO may extract wind param- 
eters from other airborne wind LIDAR systems as well. 

5 In yet another embodiment, APOLO may function by (I) 
estimating the instrument installation bias angles such as 
INS/GPS roll, pitch, and yaw, and elevation and azimuth 
biases; (2) given an airborne Doppler Wind LIDAR data and 
an INS/GPS data file, distributing the data in range gate bins; 
(3) estimating the power spectra of the monitoring pulse data 
in the beginning of each wind LIDAR data; (4) performing 
DVT (“Data Validity Test”) by investigating the peak fre- 
quency of the monitoring pulse data; (5) estimating the power 
spectra and the Doppler shift of each range gate bin data using 
their periodogram; (6) transforming the wind vector from the 
Earth coordinate (East -North-Up) into the aircraft coordinate 
using the INS/GPS data and the estimated instrument instal- 
lation bias angles; (7) estimating wind parameters in the 
20 aircraft coordinate; and (8) transforming the wind parameters 
back to the Earth coordinate. 

In a further embodiment, the wind speed and wind direc- 
tion may be computed using APOLO based on two line of 
sight (“LOS”) LIDAR returns 90 degrees apart (e.g., a line of 
25 sight at -45 and 45 degrees, where 0 degrees is straight ahead 
of the aircraft). With APOLO, simple vector summing of the 
two lines of sight may enable estimation of wind parameters. 
In this maimer, wind parameters may be estimated using 
airborne Wind LIDAR data from two orthogonal directions. 
The advantage of using only two lines of sight may be that the 
APOLO is less susceptible to noise, while the vertical wind 
component may be assumed to be zero. In another embodi- 
ment, the wind speed and wind direction may be computed 
using on a second algorithm based on five line of sight 
(“LOS”) LIDAR returns. The second algorithm may use 
pseudo-inverse to solve for unknowns (e.g., wind parameters) 
and may be sensitive to noise due to the pseudo-inverse com- 
putations. In this manner, wind parameters may be estimated 
40 using airborne Wind LIDAR data from five directions. 
Whether using two orthogonal lines of sight or five lines of 
sight, each embodiment algoritlim may compensate for 
adverse effects of different coordinate systems between the 
aircraft and the Earth. Additionally, both the embodiment two 
45 orthogonal line of sight and five line of sight algorithms may 
process wind LIDAR data in real time or near real time to 
display instantaneous wind profiles as well as long haul wind 
trend observations during an airborne mission. Additionally, 
each embodiment algorithm has the capability of instrument 
50 offset compensation, which may be critical for accurate wind 
parameter estimation. Such instrument offsets stem from a 
variety of sources such as INS/GPS unit misalignment, 
LIDAR telescope misalignment, and scanner installation 
bias. Both the APOLO algorithm and the five line of sight 
55 algorithm have a separate routine to estimate offset angles 
using ground return LIDAR data to compensate their adverse 
effects to accurate wind parameter estimates. The offset com- 
pensation routines for each algorithm may be entirely inde- 
pendent of each other except the basic principle of minimi- 
60 zation of mean square errors in a multivariate system. 

The various embodiments may address challenges faced in 
Airborne Doppler wind LIDAR systems by accounting for 
wind speeds often being far less than aircraft speeds, account- 
ing for wind directions being relative to an aircrafts naviga- 
65 tion path, accounting for INS/GPS coordinates being differ- 
ent from aircraft coordinates, accounting for instruments not 
being mounted perfectly, accounting wind definitions varying 
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by audience, accounting for difficulty in finding a reference 
wind measurement, and accounting for a lack of a reference 
algorithm. 

FIG. 1 illustrates an example aircraft mounted system suit- 
able for use with the various embodiments. The system may 5 
include an airborne Doppler Wind LIDAR system 104, an 
INS/GPS system 106, and a computing device 108 including 
a processor mounted in an aircraft 102. The aircraft 102 may 
be any type aircraft, such as a DC-8, UC- 1 2B, WB-57, Global 
Flawk Unmanned Airborne System (“UAS”), etc. The Dop- 
pier Wind LIDAR system 1 04 may be any type Doppler Wind 
LIDAR system, such as a pulsed 2-Micron Coherent Doppler 
LIDAR system. The Doppler Wind LIDAR system 104 may 
transmit andreceive pulses of light 114 along one or more line 
of sight and generate Doppler Wind LIDAR data about aero- 
sols 110 in the atmosphere. The INS/GPS system 106 may 
receive GPS timing signals from GPS satellites 112, along 
with other aircraft information, and generate INS/GPS data, 
such as the aircraft 102 Earth coordinate speed in the east, 
north, and up direction, inertial measurement unit (“IMU”) in 
meters/second, speed, roll, pitch, and heading angles in 
degrees, latitude, longitude, altitude, etc. The Doppler Wind 
LIDAR system 104 and the INS/GPS system 106 may be 
connected to the computing device 108 and may output data 
files for each line of sight (or look direction) to the computing 
device 108 which may he received by the computing device’s 
108 processor. As an example, for each line of sight along 
which data may be collected by the Doppler Wind LIDAR 
system 104, a pair of a raw Doppler Wind LIDAR data file and 
a INS/GPS data file may be output to the computing device’s 
108 processor. The Doppler Wind LIDAR data file may 
include N segments of LIDAR return data, where each seg- 
ment may be 55000 samples long resulting in a total length of 
55000’"N samples. The INS/GPS data file may also include N 
repetitions of IMU data, such as VAE, VAN, VAU, IMU speed 

(equal to ^VAE^+VAN^+VAU^), roll, pitch, and heading 
angles in degrees, longitude, latitude, and altitude. VAE, 
VAN, and VAU may be the aircraft 102 speed in meters per 
second in the Earth coordinate east, north, and up directions. 40 
At each line of sight the IMU data in the INS/GPS data file 
may be averaged over the N VAEs, VANs, VAUs, rolls, 
pitches, and heading angles and denoted as aircraft speed east 
value an aircraft speed north value (V^), an aircraft 

speed up value (V^^), an aircraft roll value (Roll), an aircraft 45 
pitch value (Pitch), and an aircraft heading value (Heading). 
The longitude and latitude at each line of sight may be the 
longitude and latitude at the first scanning or first segment. 
The Doppler Wind LIDAR data file and INS/GPS data file 
may be created and sent to the processor of the computing 50 
device 108 for use in airborne wind profiling. 

FIG. 2 illustrates a data acquisition processing and display 
system (“DAPDS”) 200 according to an embodiment. The 
DAPDS 200 may include an INS/GPS subsystem (“IGSS”) 
202 and scanner subsystem (“SCSS”) 204 connected to a 55 
processor of a computing device 210, such as core control 
subsystem (“CCSS”) or central processing unit (“CPU”) by a 
bus 205, such as a RS-232 bus. In this manner, the IGSS 202 
and/or SCSS 204 may send/receive data to/from the processor 
of the computing device 210. The DAPDS 200 may also 60 
include a digitizer subsystem (“DSS”) 206 which may receive 
analog LIDAR returns and triggers (e.g., transistor-transistor 
logic (“TTL”), etc.) and a digital signal processing (“DSP”) 
subsystem (“DSPSS”) 208. The DSS 206 and DSPSS 208 
may be connected to the processor of the computing device 65 
210 by a bus 207, such as a backplane peripheral component 
interconnect (“PCI”) bus. The processor of the computing 
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device 210 may run data acquisition and processing software 
(“DAPS”) 212 to estimate wind parameters by measuring the 
frequency shift in LIDAR back scatters while providing range 
information. 

The processor of the computing device 210 may execute 
data acquisition and processing software (“DAPS”) 212 caus- 
ing the processor of the computing device 210 to perform 
operations to use the two line of sight APOLO algorithm 
and/or a five line of sight algorithm to process wind LIDAR 
data and apply instrument offset compensation, for example 
to process the data in real time or near real time to display 
instantaneous wind profiles as well as long haul wind trend 
observations during an airborne mission. The processor of the 
computing device 210 may be connected to a keyboard/ 
mouse 218 to receive user inputs and may also be connected 
to a display 220 to display one or more plots of LIDAR data 
and/or wind parameters. The processor of the computing 
device 210 may also be connected to a memory 214 (e.g., a 
hard drive) and DVD/CD ROM 216 enabling the processor of 
the computing device 210 to read data and store data as 
needed. 

FIG. 3 illustrates example modules of data acquisition and 
processing software (“DAPS”) 301 according to an embodi- 
ment. The DAPS 301 may run within an operating system 
(“OS”) 312 nmning on a processor of a computing device. 
The DAPS 301 may include a INS/GPS control module 
(“IGCM”) 314 which may send/receive data to/from an IGSS 
302. The DAPS 301 may include a graphical user interface 
(“GUI”) control module 304 wliich may send/receive data 
to/from peripheral devices 304, such as receiving indications 
of user inputs from a keyboard or mouse and/or generating 
displays of plots and wind data for rendering on a display. The 
DAPS 302 may include a DSP control module (“DSPCM”) 
322 which may control and send/receive data to/from a 
DSPSS 310 to control a DSP. The DAPS 301 may include a 
digitizer control module (“DCM”) 320 which may send/re- 
ceive data to/from a DSS 308 to control digitization of 
LIDAR data. The DAPS 301 may include a scanner control 
module (“SCCM”) 318 which may send/receive data to/from 
a SCSS 306. The various modules 314, 316, 318, 320, 322 of 
the DAPS 301 may work in conjimction with one another to 
cause a processor of a computing device to perform opera- 
tions to use the two line of sight APOLO algorithm and/or a 
five line of sight algorithm to process wind LIDAR data with 
instrument offset compensation, for example to process the 
data in real time or near real time to display instantaneous 
wind profiles as well as long haul wind trend observations 
during an airborne mission. 

FIG. 4 illustrates an embodiment method 400 and data flow 
performed by an airborne Doppler Wind LIDAR system 
according to an embodiment. In block 402 the LIDAR system 
may generate an analog signal measurement of LIDAR 
returns. In block 404 signal conditioning techniques may be 
applied to the analog LIDAR data. In block 406 the analog 
LIDAR data may be digitized to generate a digital LIDAR 
data. In blocks 408, 410, 412, and 414 wind profiling may be 
performed using the digital LIDAR data. Specifically, in 
block 408 the digital LIDAR data may be distributed into 
range gates. In block 410 spectral analysis may be applied, for 
example to determine the power spectra and Doppler frequen- 
cies along the lines or sight. In block 412 roll, pitch, and yaw 
correction may be applied, for example to compensate for 
adverse effects of different coordinate systems between the 
aircraft and the Earth, and in block 414 a wind profiling 
algorithm, such as APOLO or a five line of sight algorithm, 
may be applied to generate a wind profile, for example by 
computing the wind speed and wind direction. In block 416 
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the wind profile may be displayed and/or archived, for 
example displayed on a display in an aircraft and/or stored in 
a memory (e.g., a hard drive) for use in post processing. 

FIGS. 5A-5D illustrate an embodiment method 500 for 
computing wind speed and wind direction from Airborne 
Wind LIDAR data. In an embodiment, the operations of 
method 500 may be performed by the processor of a comput- 
ing device in an aircraft mounted system including an air- 
borne Doppler Wind LIDAR system and an INS/GPS system. 
In block 502 of FIG. 5A the processor may receive airborne 
Doppler Wind LIDAR data. The Doppler Wind LIDAR data 
may include N segments of LIDAR returns for two or more 
lines of sight (LOS), such as two LOS perpendicular to each 
other or five LOS spaced evenly over 90 degrees. As an 
example, each segment may be 55000 samples resulting in a 
total length for the segments of the LOS being 55000*N. In 
block 504 the processor may receive INS/GPS data. The 
INS/GPS data may include N repetitions of aircraft IMU 
VAE, VAN, VAU, speed, roll, pitch, and heading angles in 
degrees. VAE, VAN, VAU may be aircraft speed in meters per 
second in the Earth coordinate east, north, and up directions. 
At each line of sight the averages of N VAEs, VANs, VAUs, 
roll, pitch, and heading angles may be denoted as aircraft 
speed east value (V^^), an aircraft speed north value (V^;,,^), an 
aircraft speed up value (V^^), an aircraft roll (Roll) value, an 
aircraft pitch (Pitch) value, and an aircraft heading value 
(Hdg). The latitude and longitude at each LOS may be the 
longitude and latitude values at the first of N scanning (i.e., 
n=0). 

In determination block 506 the processor may determine 
whether offset estimation is on or not. If offset estimation is 
on (i.e., determination block 506=“Yes”), in block 508 the 
processor may determine roll, pitch, and yaw offset angles 
based on a ground return data set. For offset estimation, in a 
manner similar to that discussed below with reference to 
block 544 of FIG. 5D the processor may find the offset angles 
in each mth range bin for m= 1 : M. One V^(x, v, t)los found, 
the processor may find Aa, A|B, A9, and AOthat minimize the 
sum of e/, where ej.=4 (Doppler shift in Hz in the kth ground 
return range bin)*L/2-IV^(x, y, L equals the 

wavelength in meters, 2.05402E-6 m and k is the index of 
ground return range bin data sets. 

In block 510 the processor may determine the number (D) 
of LOS in the LIDAR data. In block 512 of FIG. 5B the 
processor may set a value d equal to zero. In determination 
block 514 the processor may determine whether the value d is 
equal to D. In this manner, the processor may loop through 
various operations of method 500 until data from all lines of 
sight has been analyzed. If d does not equal D (i.e., determi- 
nation block 514=“No”), in block 516 the processor may 
determine the number of segments N in the LIDAR data. As 
an example, the number N may be a number of pulses aver- 
aged, such as 20. In block 518 the processor may set a value 
of n equal to zero. In determination block 520 the processor 
may determine whether n equals N. In this manner, the pro- 
cessor may loop through various operations of method 500 
until data from all segments of the LOS has been analyzed. 

Ifndoesnot equalN(i.e., determination block 520=“No”), 
in block 522 of FIG. 5C the processor may sort the nth 
segment into a number of range gate bins M. As an example, 
the processor may acquire s„(t) and divide it into M+1 vec- 
tors, Xq^„, Xj „ . . . X^„, with an overlap excluding the first 
vector Xq „. (Xq „ may not have an overlap with X j „ and may 
only be used for a data validation test as discussed below). 
Note that the length of the first vector Xq „ may be different 
from the rest of the vectors, for example, the first vector may 
be 1024 samples long while the rest have 512 samples. The 


8 

first vector may include pre-triggering data (1,000 samples 
for example were used in the Genesis and Rapid Intensifica- 
tion Processes (“GRIP”) test mission discussed below) as 
well. In such a case, the total number of data, such as 55,000, 
5 may include the pre-triggering data. In block 524 the proces- 
sor may determine the power spectra of each range gate bin 
according to the equation X„, „=Power spectrum of x„, Vm. 

In block 526 the processor may perform a data validity test 
to find a reference signal and count the number of passing 
to pulses (Q). As an example, the data validity test (“DVT”) may 
be to check if ZDmin<fQ „<ZDmax, where 


hn =max[Xo^(/)]. 

15 

Note that maximum power frequency may be estimated 
above a nominal value F Hertz (“Hz”). ZDmin and ZD max 
may be the minimum and maximum limits of Doppler fre- 
20 quency, which were 95 MHz and 115 MHz, respectively, 
during the GRIP test missions discussed below. If true and it 
is the first case, will be the reference zero Doppler fre- 
quency. If the current fg^ does not pass the DVT, go to the next 
(n+l)th interation and test it. Assuming the fg^ passes the 
25 DVT for the first time at the Kth iteration, where n=K which 
can be anywhere between 0 and N- 1 , then s^(t) becomes the 
reference signal and is the reference zero Doppler fre- 
quency for all n at the current look direction d. In block 528 
the processor may determine the frequency change Afg „ as 
30 equal to fo^-fo,„, naK. 

In block 530 the processor may determine the zero Doppler 
alignment (“ZDA”) X„,„(f) as equal to X„,„(f+Afo,„), for 
m=l:M. Note that m does not start at 0, meaning the first 
vector Xo_„ may be excluded in the data processing since Xq,„ 
35 is used only to estimate fo M is the number of range bins. 
Also, note that the zero Doppler alignment in the airborne 
system differs from that of ground based wind profiling algo- 
rithms which use X„_„(f-Afo„). For example, assume 
fo,*: =100 MHz from Xg j; in s^(t) and fo^^=98 MHz from Xq,^ in 
40 s^(t), where K<q<N - 1 , and the LIDAR returned data at n=K, 
s^t), passed the DVT for the first time, and the later pulse 
return at n=q, s^(t) also passed the DVT. Then Afo^^=fo^- 
ig^=2 MHz resulting inX„_^(f)=X„,^^(f+2 MHz), form=l :M. 
This means shifting X„ ,^(1) by 2 MHz to the left (not the right) 
45 to obtain X„_^(f). This differs from the ground based wind 
profiling algoritlun where the shift is to the right. In block 532 
the processor may increment n. 

In FIG. 5B, when n equals N (i.e., determination block 
520=“Yes”). in FIG. 5D bock 534 the processor may deter- 
50 mine the power frequency for each of the data segment as a 
value Y equal to { Yj . . . Y^} wherein 


1 0 “' 

^n=0 

Q is the number of pulse returns having passed the DVT 
(assuming Q consecutive raw data arrays have passed the 
00 p)VT). Inblock536theprocessormay determine the Doppler 
shift f as a value equal to {fj • • • ^m} wherein 


fm= foK + max[KOT(/)] - 100 MHz, m = 1 : M. 
65 ’ f>F 
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Note that the maximum power frequency is estimated above 
a nominal value F Flz and includes a sign change and 100 
MFlz term when compared with the ground based wind pro- 
filing algorithm 


fm = fo.K -max[y„(/)]. 
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6B.) In each mth range bin the following may be then be 
computed for m=l :M. First v^(N,E:U) may be defined as: 


5 

Van 


^| '^AE + '^AN 

Vy,(/V, E, U) = 

VaE 

= 

0 


Vau . 




The 100 MFIz term may be of a different value for different 
LIDAR instrument designs. In block 538 the processor may 
determine the power spectra p as a value equal to {pj . . . p^,^} 
wherein p„=max[Y„(f)], f>F, and m=l :M. 

In block 540 the processor may determine the LOS range as 
a value equal to Bin Range [meter] *(1 -Overlap Sample 
Count/Bin Sample Count)*m+Bin Range [meter]*0.5+ 
r*c*0.5/ F^ where m=0:M-l, Overlap Sample Count is the 
number of samples in the overlap between adjacent range 
bins, r is the index of the first sample in Xj c is the speed of 
light in meters per second, F^ is the sample frequency in Hz, 
and Bin Range[meter]=Bin Sample Count*c*0.5/F„. (The 
distance between each sample is c*0.5/F^.) Note r in the last 
term must take into consideration the length of pre-triggering 
data. For example, if the first 1024 samples were used for the 
DVT and the length of the pre-triggering data was 1000, r 
would be 24 (i.e., 1024-1000) instead of 1024. Note that the 
index begins at zero. If the length of pre-trigger was zero (i.e., 
no pre-triggering data), then r would be 1024. The last term 
r*c*0.5/Fj is added to apply the offset caused by excluding 
the first bin data less the pre-triggering data in data process- 
ing. The second term Bin Range [meter] *0.5 is added in order 
to go by the center of each bin range rather than the beginning 
or the end. Since Xq „ is used only for the DVT, the second bin 
data Xi „ must not include samples from Also note that 
the length of the first and the second bin data will likely be 
different such as 1 024 for the first bin and 500 or 5 1 2 for the 
rest. Therefore, the data index r of the first sample in the 
second bin x^ „ must not be one of the indices of the data in 
Xq,„. Again, the value of r must be adjusted if pre-trigger data 
are used for the DVT. In block 542 the processor may plot the 
LOS range vs f and p (i.e., the line of sight range versus the 
Doppler Shift and Power spectra). The LOS range may be 
zero at the aircraft and increases as the LIDAR beam travels 
away from the aircraft. 

In block 544 the processor may determine speed east 
(V^), speed north (V wiv)’ speed up (ypyu) wind com- 
ponents for each line of sight (d) as vectors defined by: 


VwiN, E, = 


The wind components speed east (V^^), speed north 
(Vrw): and speed up (Vy^j) are determined by converting tan 



to an angle from N and determining an aircraft yaw (Yaw) 
value as the aircraft heading value minus the angle from N. 
(An illustration of such relationships is illustrated in FIG. 


Note that v^, v^^, and v^^may be on a “flat” plane since 
R(Roll), P(Pitch), and Y (Yaw) matrices have not been applied 
yet. Also, note that v^^^ is now aligned with the nose of the 
aircraft. Second, V^(N, E, U) may be found as: 


VA(x,y,z) = 




= R(a + Att)P(/? + A/?)y(7 + Ay)V.,(A', E,-U) = 


20 


R(a + ^a)P(P + ^P)Y(y + Ay) 


V + ^AN 


0 

-Vau 


25 where 


R{a) = 


10 0 ■ 
0 cosa sino' , 
0 -sina cosa 


30 


cos/9 

0 

-sin/9' 


cosy 

siny 

O' 


p{p) = 

0 

1 

0 

“r: 

II 

-siny 

cosy 

0 



sin/9 

0 

cos/9 


0 

0 

1 . 


35 These matrices rotate the axes of the aircraft while the 
vector remains in the same positions. Aa, Ap, and Ay are the 
roll, pitch, and yaw offset values that need to be pre-calcu- 
lated. Now x, y, and z axes (and the aircraft) are crooked since 
roll, pitch, and yaw have been applied. The vector v^(N, E, 
40 -U) remains in the same position, but the aircraft axes (x, y, 
and z, illustrated for general reference to an aircraft in FIG. 
6A and FIG. 6C illustrates the relationships of pitch, roll, and 
yaw to aircraft orientation) have been rotated. The values 

and V^ are in reference to the new (crooked) x, y, and z 
45 coordinates. Third, V^(x, y, z)^q^ may be found which is 
V^(x, y, z) along the LOS by projecting V^(x, y, z) onto the 
unit vector Ugjj^x, y, z) along the LOS(0+A0, A0), where 

0=6O°-O.12°(offset)=59.88° (Elevation: x to z) and 0=-45°, 
-22.5°, 0°, 22.5°, 45° (d=0:4) (Azimuth: x to y) A0 and A0are 
50 the elevation and the azimuth offset relative to the aircraft 
body. V^(x, y, z)^os^(''. 4 'Ue o)Ue » where • indicates inner 
product, and 


55 


“0.0 (J:, y, 4) = “0.0 


cos(0 + A^)cos(0 + A0) 


cos(0 + A0)sin(0 + A0) 


= LOS unit vector. 


sin(0 + A0) 


60 Fourth, the wind magnitude IV ^^x, y, z)l along the LOS(0, 
0 ) may be found such that: IV^x, y, z)^o^l=IV^x, y, z)+V^ 
(x, y, z)l along the LOS-IV^(x, y, zj^^^Nf^ (Doppler shift in 
Hz in the mth range bib)*L/2-IV^(x, y, where L 

equals the wavelength in meters, 2.05402E-6 m. Fifth, the 
65 IV^x,y,z),^p;^„gl=IV^x,y,z)l in the xy plane of the aircraft 
(which is not level) may be found. This is not the projection 
onto the xy plane since the value before being projected is 
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12 


desired, which is larger than IV^x , y, z)l. Thus, IV^x, y, 
z)xj>^>za«el=IV^x, y, z)l along the LOS(6, 0) discussed above 
in the fourth step 


1 

cos(0 + A0) 


= \Vw(x, y, z)ios\'^ 


1 

cos(^ + A0)' 


Sixth, IV^x, y, may be decomposed into its com- 

ponents in the x, y, and z directions to construct the wind 
vector V jf,(x, y, z). Assuming the z component is zero (i.e., 
that there is no vertical wind component at the aircraft in the 
crooked position). 





X cos(0 + A0) 

Vw{x, y, z) = 

'^Wy 

= 

\Vwix, y, xsin(0 + A0) 




0 


compute the wind components using two lines of sight per- 
pendicular to each other. In another embodiment, the proces- 
sor may apply a five line of sight based algorithm that uses 
five lines of sight equally spaced over a 90 degree arc to 
5 compute the wind speed and the wind direction using pseudo- 
inverse functions. 

The following illustrates an embodiment method for find- 
ing wind speed and wind direction using the APOLO. In each 
range bin (m), the north, east, and up wind components may 
be found using E, at the first and the fifth LOSs 

(e.g., d=l and 5) as follows: 


15 


Vw(N, E, U)„ = 


VwEjn 

VwU,m 




\yidl) ] 
'^WN.m 


r 1 

^WN.m 


vw(N, E, Ufy + vw(N, E. Ufy = 

w(l) 

^WE.m 

+ 

w(5) 

^WE.m 

20 


w(l) 

L '^WUjn J 


w(5) 

L '^WU.m J 


Seventh, 


Vw(N,E, U) = 


^WN 

VWE 


^WU 


may be found from V,p(x, y, z) by computing 


Then, the horizontal wind speed (“HWS”) in the mth range 
bin may be obtained according to: 

25 

t'w/V.m ^WE.m = 'J /n + ^WN ,m) + (^WE.m + ^WE.m) ’ 

m = l:M 
30 


Vw(N, E, -U) = 


^WN 

y WE 


V wu 


[/?(ar + Aa)P(^ + + Ay)] ‘ Vw(x, y, z)- 


To compute the horizontal wind direction (“HWD”), first 
may be computed as follows: 


35 



Here, V^, V^, and are back on a “flat” or leveled xy 
plane of the aircraft without roll, pitch, or yaw, and the x axis 
(nose of the aircraft) aligns with the direction of as 
illustrated in FIG. 6D. may equal tan 



and 72 may equal yy 2=7 1 - 90° . Then, V ^N, E, U) (or v E, 

U)„^‘^ in the mth range bin at the dth LOS may be found by: 



VWN 


Vwn^mji + Vvf£siny2 

Vw(A', E, U) = 

'^WE 

= 

Viv^cosyi +Vi4/£COS72 


. '^WU . 


-{-Vwu) 


Vw(N, E, Uff = 




Not that 8„ is measured counter-clockwise as positive 
starting positive east (positive east=0°, positive north=90°, 
negative east=180°, negative north=270°, etc.). The math- 
ematical HWD in the mth range bin is determined as follows, 
if 8„<0, (HWD„)„,,,=360°+6„, else (HDW„)„,,,,,„ , where 
m=l:M. The meteorological HWD in the mth range bin is 
45 determined as follows, if 90°-(HWD„,)„,^^;,<0, 
HWD„=360°+90°-(HWD„)„^,;,, else HWD„=90°- 
(HWD„,)„,^;^, where m=l :M. Or simply HWD„,=mod(90°- 
(HWD„)„^,;,),360°) where mod (x,y)=x-y*floor(x/y)=re- 
mainder of x/y. 

50 In block 550 altitude may be computed and stored. As an 
example altitude may be computed using 0=6O°-O.12°(off- 
set)=59.88°. Altitude may equal GPS aircraft altitude 
[meters]-LOS Range [meter] *cos(9O°-0) which may equal 
Altitude [meter] -LOS Range [meter] *cos(30. 12°). LOS 
55 range may have been previously calculated as discussed 
above, and an assumption that there is zero roll, pitch, and 
yaw angles may be made, otherwise LOS range may not be 
converted uniformly into altitudes. In block 552 the processor 
may plot and store wind speed and wind direction versus 
60 altitude, for example plotting HWS vs Altitude and HWD vs 
altitude on one or more displays. 


In block 546 the processor may increment d. If d equals D 
(i.e., determination block 514 in FIG. 5B=“Yes”), in block 
548 the processor may compute and store wind speed and 65 
wind direction using the north, east and up wind components. 

In an embodiment, the processor may apply the APOLO to 


Example Results of Wind Profiling Using the 
Various Embodiments 

The various embodiments have been tested in fiights using 
a pulsed 2-micron coherent Doppler LIDAR system and a 
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data acquisition, processing, and display system (“DAPDS”) 
which estimates wind parameters by measuring the frequency 
shift in LIDAR back scatters while providing range informa- 
tion. A test system flew during the GRIP hurricane campaign 
on a DC-8, which flew out of Fort Lauderdale, Fla. The total 
flight time was approximately 130 hours. The goal was to test 
the airborne wind profiling capabilities of the system. A sec- 
ond flight demonstration and a third flight demonstration on a 
UC-12B were also conducted. Flight tests on the UC-12B 
centered on evaluating the installation of the system into a 
smaller aircraft and in the high vibration environment asso- 
ciated with turboprop engines . The system was operated in 22 
and 25 hours of flights on the UC-12B, respectively. The 
following discussion of FIGS. 7-15 discuss example results 
gathered in the various test flights using the various embodi- 
ments. 

FIG. 7 shows the result of offset compensation in horizon- 
tal wind speed (“FIWS”) from APOLO on a DCS while it was 
at 10,608 meters in altitude before (on the left) and after (on 
the right) applying offset compensation. The APOLO based 
estimates are shown as the circles in the graphs. A total of 5 1 2 
samples were used to represent wind at each altitude, and a 
4,096-point Fast Fourier Transform (“FFT”) was used to pro- 
cess each set of 512 samples using zero padding. It clearly 
shows the effectiveness of offset compensation with APOLO. 
The APOLO estimates near 4.5-9 km in altitude are erroneous 
due to low received signal backscatter associated with atmo- 
spheric conditions of low aerosol content. The solid line is the 
dropsonde measurement only for a reference purpose, and not 
to imply which measurements are more accurate (the two 
sensors measure different volumes over different time 
scales). The good agreement on the right plot means that the 
wind was slowly varying over the dropsonde fall time and 
overthe distance of DC-8 movement for the LIDAR measure- 
ment time. This agreement is not guaranteed in a general wind 
field. When it does fortunately occur, it permits comparison of 
processing algorithms. FIG. 8 shows the result of offset com- 
pensation in horizontal wind direction (HWD) from APOLO. 
The data are the same as FIG. 7. The adverse effects of a low 
Signal-to-Noise-Ratio (SNR) are clear near 4.5-9 km in alti- 
tude. FIG. 9 shows the result of offset compensation in hori- 
zontal wind speed (“HWS”) from the five line of sight algo- 
rithm before (on the left) and after (on the right) applying 
offset compensation. The five line of sight algorithm based 
estimates are shown as the circles in the graphs. A total of 5 1 2 
samples were used to represent wind at each altitude, and a 
4,096-point Fast Fourier Transform (“FFT”) was used to pro- 
cess each set. The solid line shows dropsonde measurements 
for reference purposes. FIG. 10 shows the result of offset 
compensation in horizontal wind direction (HWD) from the 
five line of sight algorithm. The data are the same as FIG. 9. 

However, HWS and HWD of the second five line of sight 
algorithm are compared with those of APOLO in FIGS. 11 
and 12. FIG. 11 illustrates the HWS from APOLO in the left 
and the second five line of sight algorithm on the right. Offset 
compensation was applied in both plots. The data was taken 
during the GRIP mission ontheDC8. FFT size was 4096 with 
zero padding. Altitude=I 0,608 meters. The solid line repre- 
sents dropsonde data, and dots are the estimates from APOLO 
(left) and the second five line of sight algorithm (right). FIG. 
11 compares HWS from APOLO (left) with that from the 
second five line of sight algorithm (right). The same data were 
used in this comparison, and dropsonde measurements are 
shown for a reference purpose again. FIG. 12 compares HWD 
in the same fashion. FIG. 12 shows HWD from APOLO in the 
left graph and the second five line of sight algorithm in the 
right graph. Offset compensation was applied during the 
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GRIP mission on the DC8. FFT size was 4096 with zero 
padding. Altitude was 10,608 meters. The line represents 
dropsonde data, and black dots are the estimates from 
APOLO on the left and the second five line of sight algorithm 
5 on the right. 

FIG. 13-15 show panoramic displays of power, HWS, and 
HWD during the flight demonstration in a UC-12B aircraft 
and flew at lower altitudes than the GRIP hurricane campaign. 
The time stamps are in Zulu time and the images were gen- 
10 erated in software developed for post processing to create a 
long-haul observation displays as shown. FIG. 13 is a pan- 
oramic display of power distribution from the flight campaign 
on the UC-12B. The algorithm used was APOLO with offset 
compensation applied and a FFT size of 512. FIG. 14 is a 
15 panoramic display of HWS from the flight campaign on the 
UC-12B. The algorithm used was APOLO with offset com- 
pensation applied and a FFT size of 512. FIG. 15 is a pan- 
oramic display of HWD from the flight campaign on the 
UC-12B. The algorithm used was APOLO with offset com- 
20 pensation applied and a FFT size of 512. The aircraft flight 
track on the day the data in FIGS. 13-15 was gathered was 
from Hampton, Va. to Charlotte, N.C. and back. While over 
Charlotte the aircraft made multiple oveipasses of a ground- 
based wind profiling radar for comparing the measurements 
25 of the two sensors. The power spectral density time trend of 
FIG. 13 reveals cloud and aerosol structure through the course 
of the flight with an aerosol layer evident from 1 .5 to 3.5-km 
altimde throughout most of the flight. The horizontal wind 
speed trend of FIG. 14 shows a horizontal wind magnitude 
30 vertical shear at around 1 -km altitude, typical of the transition 
from the atmospheric boundary layer (where the wind speed 
is about 5-m/s in this case) to the free troposphere (where the 
wind speed is about 1 5-m/s in this case). The wind direction 
plot of FIG. 15 shows a wind direction from the northwest in 
35 the free troposphere above 1 -km altitude, but a more south- 
erly direction in the low altitudes in the more turbulent atmo- 
spheric boundary layer below 1-km altitude, a horizontal 
wind direction vertical shear. 

The various embodiments may be implemented in any of a 
40 variety of computer devices, an example of which is illus- 
trated in FIG. 16. A computer device 1600 will typically 
include a processor 1601 coupled to volatile memory 1602 
and a large capacity nonvolatile memory, such as a disk drive 

1605 of Flash memory. The computer device 1600 may also 
45 include a floppy disc drive 1603 and a compact disc (CD) 

drive 1604 coupled to the processor 1601. The computer 
device 1600 may also include a number of connector ports 

1606 coupled to the processor 1601 for establishing data 
connections or receiving external memory devices, such as a 

50 USB or FireWire® connector sockets, or other network con- 
nection circuits for coupling the processor 1601 to a network 
or bus. The computer device 1600 may also include the track- 
ball 1607, keyboard 1608 and display 1609 all coupled to the 
processor 1601. 

55 The processor 1601 may be any programmable micropro- 
cessor, microcomputer or multiple processor chip or chips 
that can be configured by software instructions (applications) 
to perform a variety of functions, including the functions of 
the various embodiments described above. In some devices, 
60 multiple processors may be provided, such as one processor 
dedicated to wireless communication functions and one pro- 
cessor dedicated to running other applications. Typically, 
software applications may be stored in the internal memory 
1602 and/or 1605 before they are accessed and loaded into the 
65 processor 1601. The processor 1601 may include internal 
memory sufficient to store the application software instruc- 
tions. In many devices the internal memory may be a volatile 
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or nonvolatile memory, such as flash memory, or a mixture of 
both. For the purposes of this description, a general reference 
to memory refers to memory accessible by the processor 1601 
including internal memory or removable memory plugged 
into the device and memory within the processor 1601 itself. 

The foregoing method descriptions and the process flow 
diagrams are provided merely as illustrative examples and are 
not intended to require or imply that the steps of the various 
embodiments must be performed in the order presented. As 
will be appreciated by one of skill in the art the order of steps 
in the foregoing embodiments may be performed in any order. 
Words such as “thereafter,” “then,” “next,” etc. are not 
intended to limit the order of the steps ; these words are simply 
used to guide the reader through the description of the meth- 
ods. Further, any reference to claim elements in the singular, 
for example, using the articles “a,” “an” or “the” is not to be 
construed as limiting the element to the singular. 

The various illustrative logical blocks, modules, circuits, 
and algorithm steps described in connection with the embodi- 
ments disclosed herein may be implemented as electronic 
hardware, computer software, or combinations of both. To 
clearly illustrate this interchangeability of hardware and soft- 
ware, various illustrative components, blocks, modules, cir- 
cuits, and steps have been described above generally in terms 
of their functionality. Whether such functionality is imple- 
mented as hardware or software depends upon the particular 
application and design constraints imposed on the overall 
system. Skilled artisans may implement the described func- 
tionality in varying ways for each particular application, but 
such implementation decisions should not be interpreted as 
causing a departure from the scope of the present invention. 

The hardware used to implement the various illustrative 
logics, logical blocks, modules, and circuits described in 
connection with the aspects disclosed herein may be imple- 
mented or performed with a general purpose processor, a 
digital signal processor (DSP), an application specific inte- 
grated circuit (ASIC), a field programmable gate array 
(FPGA) or other programmable logic device, discrete gate or 
transistor logic, discrete hardware components, or any com- 
bination thereof designed to perform the functions described 
herein. A general-purpose processor may be a microproces- 
sor, but, in the alternative, the processor may be any conven- 
tional processor, controller, microcontroller, or state 
machine. A processor may also be implemented as a combi- 
nation of computing devices, e.g., a combination of a DSP 
and a microprocessor, a plurality of microprocessors, one or 
more microprocessors in conjunction with a DSP core, or any 
other such configuration. Alternatively, some steps or meth- 
ods may be performed by circuitry that is specific to a given 
function. 

In one or more exemplary aspects, the functions described 
may be implemented in hardware, software, firmware, or any 
combination thereof. If implemented in software, the func- 
tions may be stored as one or more instructions or code on a 
non-transitory computer-readable medium or non-transitory 
processor-readable medium. The steps of a method or algo- 
rithm disclosed herein may be embodied in a processor-ex- 
ecutable software module which may reside on a non-transi- 
tory computer-readable or processor-readable storage 
medium. Non-transitory computer-readable or processor- 
readable storage media may be any storage media that may be 
accessed by a computer or a processor. By way of example but 
not limitation, non-transitory computer-readable or proces- 
sor-readable media may include RAM, ROM, EEPROM, 
FLASH memory, CD-ROM or other optical disk storage, 
magnetic disk storage or other magnetic storage devices, or 
any other medium that may be used to store desired program 


16 

code in the form of instructions or data structures and that 
may be accessed by a computer. Disk and disc, as used herein, 
includes compact disc (CD), laser disc, optical disc, digital 
versatile disc (DVD), floppy disk, and blu-ray disc where 
5 disks usually reproduce data magnetically, while discs repro- 
duce data optically with lasers . Combinations of the above are 
also included within the scope of non-transitory computer- 
readable and processor-readable media. Additionally, the 
operations of a method or algorithm may reside as one or any 
to combination or set of codes and/or instructions on a non- 
transitory processor-readable medium and/or computer-read- 
able medium, which may be incorporated into a computer 
program product. 

The preceding description of the disclosed embodiments is 
15 provided to enable any person skilled in the art to make or use 
the present invention. Various modifications to these embodi- 
ments will be readily apparent to those skilled in the art, and 
the generic principles defined herein may be applied to other 
embodiments without departing from the spirit or scope of the 
20 invention. Thus, the present invention is not intended to be 
limited to the embodiments shown herein but is to be 
accorded the widest scope consistent with the following 
claims and the principles and novel features disclosed herein. 

25 What is claimed is: 

1. An airborne wind profiling method performed by an 
aircraft mounted system comprising an airborne Doppler 
Wind LIDAR system, a INS/GPS system, and a computing 
device connected to the airborne Doppler Wind LIDAR sys- 
30 tern and the INS/GPS system, the computing device compris- 
ing a display, a memory, and a processor connected to the 
display and the memory, the method comprising: 

receiving airborne Doppler Wind LIDAR data and associ- 
ated INS/GPS data, wherein the Doppler Wind LIDAR 
35 data includes data segments for two or more lines of 
sight; 

distributing the data segments of each line of sight into 
range gate bins; 

determining a power spectra of each range gate bin; 

40 performing a data validity test based at least in part on the 
power spectra of each range gate bin to determine a 
reference signal for each of the data segments; 

determining a frequency change from the reference signal 
for each of the data segments; 

45 determining a zero Doppler alignment for each of the data 

segments based on the frequency change; 

determining a power frequency for each of the data seg- 
ments based on the data segment’s zero Doppler align- 
ment; 

50 determining a Doppler shift f for each line of sight based on 
the line of sight’s respective maximum power fre- 
quency; 

determining a power spectra p for each line of sight based 
on the line of sight’s respective maximum power mag- 
55 nitude; 

determining a line of sight range for each line of sight; 

plotting the line of sight range for each line of sight versus 
the line of sight’s Doppler shift f and power spectra p; 

determining east, north, and up wind components for each 
60 line of sight based on the plot of the line of sight range for 

each line of sight versus the line of sight’s Doppler shift 
f and power spectra p; 

computing a wind speed and a wind direction using the 
east, north, and up wind components for each line of 
65 sight; 

computing an altitude based on the INS/GPS data and line 
of sight range for each line of sight; and 
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plotting the wind speed versus the altitude and the wind 
direction versus the altitude on the display. 

2. The method of claim 1, wherein determining a power 
frequency for each of the data segments based on the data 
segment’s zero Doppler alignment comprises: 

counting a number of pulse returns Q that have passed the 
data validity test; and 

determining the power frequency for each of the data seg- 
ments as a value Y equal to { Yj . . wherein 


1 


M is a number of range gate bins, n is the data segment 
number, and „(f) is the zero Doppler alignment. 

3. The method of claim 2, wherein: 
determining a Doppler shift f for each line of sight based on 
the line of sight’s respective maximum power frequency 
comprises determining the Doppler shift f as a value 
equal to {fi . . . wherein 


= h.K +max[t'„(/)] - 

fzF 


100 MHz, m = \:M, 
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fo,^: is the reference signal, and F is a nominal value in 
Hz; and 30 

determining a power spectra p for each line of sight based 
on the line of sight’s respective maximum power fre- 
quency comprises determining the power spectra p as a 
value equal to {p^ . . . pJ^^} wherein 
;j„=max[7„(/)],/£F, andm = l:M. 35 

4. The method of claim 3, wherein the associated INS/GPS 
data includes an aircraft speed east value an aircraft 

speed north value an aircraft speed up value (V^u), an 

aircraft roll value, an aircraft pitch value, and an aircraft 
heading value for each segment for the two or more lines of 40 
sight, and 

wherein determining east, north, and up wind components 
for each line of sight based on the plot of the line of sight 
range for each line of sight versus the line of sight’s 
Doppler shift f and power spectra p comprises determin- 45 
ing speed east (V^^), speed north and speed up 

(Swu) wind components for each line of sight (d) as 
vectors defined by: 


translating the aircraft vector for each range gate bin 
in reference to an adjusted aircraft axis using the 
aircraft speed east value (V^^), the aircraft speed 
north value (V^), the aircraft speed up value 
the aircraft roll value, the aircraft pitch 
value, the aircraft heading value, the aircraft yaw 
value, and roll, pitch, yaw, elevation, and azimuth 
offset values; 

estimating wind (W) components in reference to the 
adjusted aircraft axis based on the Doppler shift f 
for each range gate bin and the translated aircraft 
vector for each range gate bin; and 
translating the wind (W) components in reference to 
the adjusted aircraft axis into Earth reference wind 
components speed east (S we)^ speed north (Vwn), 
and speed up (yfrv)- 

5. The method of claim 4, wherein the roll, pitch, yaw, 
elevation, and azimuth offset values are estimated at least in 
part on a Doppler shift in a reference ground return range bin 
and a translated aircraft vector for the reference ground return 
range bin. 

6. Tlie method of claim 5, wherein the Doppler Wind 
LIDAR data includes data segments for only two lines of sight 
which are perpendicular to each other, and wherein comput- 
ing a wind speed and a wind direction using the east, north, 
and up wind components for each line of sight comprises: 

adding the wind components for each line of sight such 
that: 

the wind speed is equal to: 





and 

the wind direction is equal to: 




where 




= 360° + (5„ if (5„ < 0 othenvise = S„ 


-2 = 1 : M and = tan 


-1 ''WN,n. 

1/(1) 

V ^WE,rr 




+ V'w 


vw(/V, E, = 


50 


55 


Wherein wind components speed east (V^), speed 
north (y wn)^ ond speed up (V yyjj) are determined by: 
converting 

60 


to an angle from N; 65 

determining an aircraft yaw value as the aircraft head- 
ing value minus the angle from N; 


7. An aircraft mounted system, comprising: 
an airborne Doppler Wind LIDAR system; 

INS/GPS system; and 

a computing device connected to the airborne Doppler 
Wind LIDAR system and INS/GPS system, the comput- 
ing device comprising: 
a display; 
a memory; and 

a processor connected to the display and the memory, 
wherein the processor is configured with processor- 
executable instructions to perform operations com- 
prising: 

receiving airborne Doppler Wind LIDAR data from 
the airborne Doppler Wind LIDAR system and 
associated INS/GPS data from the INS/GPS sys- 
tem, wherein the Doppler Wind LIDAR data 
includes data segments for two or more lines of 
sight; 
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distributing the data segments of each line of sight 
into range gate bins; 

determining a power spectra of each range gate bin; 
performing a data validity test based at least in part on 
the power spectra of each range gate bin to deter- 
mine a reference signal for each of the data seg- 
ments; 

determining a frequency change from the reference 
signal for each of the data segments; 
determining a zero Doppler alignment for each of the 
data segments based on the frequency change; 
determining a power frequency for each of the data 
segments based on the data segment’s zero Doppler 
alignment; 

determining a Doppler Shift f for each line of sigjit 
based on the line of sight’s respective maximum 
power frequency; 

determining a power spectra p for each line of sight 
based on the line of sight’s respective maximum 
power magnitude; 

determining a line of sight range for each line of sight; 
plotting the line of sight range for each line of sight 
versus the line of sight’s Doppler shift f and power 
spectra p on the display; 

determining east, north, and up wind components for 
each line of sight based on the plot of the line of 
sight range for each line of sight versus the line of 
sight’s Doppler shift f and power spectra p; 
computing a wind speed and a wind direction using 
the east, north, and up wind components for each 
line of sight; 

computing an altitude based on the INS/GPS data and 
line of sight range for each line of sight; and 
plotting the wind speed versus the altitude and the 
wind direction versus the altitude on the display. 

8. The system of claim 7, wherein the processor is config- 
ured with processor executable instructions to perform opera- 
tions such that determining a power frequency for each of the 
data segments based on the data segment’s zero Doppler 
alignment comprises: 

counting a number of pulse returns Q that have passed the 
data validity test; and 

determining the power frequency for each of the data seg- 
ments as a value Y equal to { Yj . . . Y^,^} wherein 


1 


m = 1 : M, 


M is a number of range gate bins, n is the data segment 
number, and X„ „(f) is the zero Doppler alignment, 

9. The system of claim 8, wherein the processor is config- 
ured with processor-executable instructions to perform 
operations such that: 

determining a Doppler shift f for each line of sight based on 
the line of sight’s respective maximum power frequency 
comprises determining the Doppler shift f as a value 
equal to {fj • • • u wherein 


fm = k.K +max[y„(/)] - 

f-sF 


100 MHz, m = 1 : M, 


is the reference signal, and F is a nominal value in 
Hz; and 

determining a power spectra p for each line of sight based 
on the line of sight’s respective maximum power fre- 
5 quency comprises determining the power spectra p as a 
value equal to {p^ . . . p^^} wherein 

/)„=max[}^(/)],^F, and m=l \M. 

10. The system of claim 9, wherein the processor is con- 
ic figured with processor-executable instructions to perform 
operations such that: 

the associated INS/GPS data includes an aircraft speed east 
value (V^^), an aircraft speed north value (V^jv), an 
aircraft speed up value an aircraft roll value, an 

15 aircraft pitch value, and an aircraft heading value for 
each segment for the two or more lines of sight, and 

determining east, north, and up wind components for each 
line of sight based on the plot of the line of sight range for 
each line of sight versus the line of sight’s Doppler shift 
20 f and power spectra p comprises determining speed east 
(Vwe\ speed north (V^), and speed up (V^) wind 
components for each line of sight (d) as vectors defined 
by: 


25 


30 


35 


vw{N,E, = 


'(d) ■ 

WNjn 


id) 

WEjn 


id) 

WUjn 


wherein wind components speed east (V^), speed 
north (Vb 7 v,)j stid speed up (Vjyu) sr® determined by 
converting 


C 


40 to an angle from N; 

determining an aircraft yaw value as the aircraft head- 
ing value minus the angle from N; 
translating the aircraft vector for each range gate bin 
in reference to an adjusted aircraft axis using the 
45 aircraft speed east value (V^^), the aircraft speed 

north value (V^), the aircraft speed up value 
the aircraft roll value, the aircraft pitch 
value, the aircraft heading value, the aircraft yaw 
value, and roll, pitch, yaw, elevation, and azimuth 
50 offset values; 

estimating wind (W) components in reference to the 
adjusted aircraft axis based on the Doppler shift f 
for each range gate bin and the translated aircraft 
vector for each range gate bin; and 
55 translating the wind (W) components in reference to 

the adjusted aircraft axis into Earth reference wind 
components speed east (V^), speed north (V^), 
and speed up (Vfyu)- 

11. The system of claim 10, wherein the processor is con- 
60 figured with processor-executable instructions to perform 
operations such that the roll, pitch, yaw, elevation, and azi- 
muth offset values are estimated at least in part on a Doppler 
shift in a reference ground return range bin and a translated 
aircraft vector for the reference ground return range bin. 

65 12. The system of claim 11, wherein the processor is con- 

figured with processor-executable instructions to perform 
operations such that the Doppler Wind LIDAR data includes 



US 9,007,570 B1 


21 

data segments for only two lines of sight which are perpen- 
dicular to each other, and wherein computing a wind speed 
and a wind direction using the east, north, and up wind com- 
ponents for each line of sight comprises: 

adding the wind components for each line of sight such 5 
that: 

the wind speed is equal to: 

V (C,™ + + (vSL + v^weJ . 

m = 1 : M\ 

and 15 

the wind direction is equal to: 
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computing an altitude based on the INS/GPS data and line 
of sight range for each line of sight; and 

plotting the wind speed versus the altitude and the wind 
direction versus the altitude on a display connected to 
the processor. 

14. The non-transitory processor readable medium of 
claim 13, wherein the stored processor-executable instruc- 
tions are configured to cause the processor to perform opera- 
tions such that determining a power frequency for each of the 
data segments based on the data segment’s zero Doppler 
alignment comprises: 

counting a number of pulse returns Q that have passed the 
data validity test; and 

determining the power frequency for each of the data seg- 
ments as a value Y equal to {Yj . . . Yj^} wherein 


mod(90° - 360”), 

where 

{HWDm)moth = 360° + (5^ if (5^ < 0 otherwise 


m = l\M and 6^ 


- tan 


^w/v,» 

+ '^WN,m 

V 

, w(5) • 

+ ''WE,m / 
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13. A non-transitoiy processor readable medium having 
stored thereon processor-executable instructions configured 
to cause a processor of a computing device of an aircraft 
mounted system comprising an airborne Doppler Wind 
LIDAR system and a INS/GPS system connected to comput- 
ing device to perform operations comprising: 

receiving airborne Doppler Wind LIDAR data and associ- 
ated INS/GPS data, wherein the Doppler Wind LIDAR 35 
data includes data segments for two or more lines of 
sight; 

distributing the data segments of each line of sight into 
range gate bins; 

determining a power spectra of each range gate bin; 40 
performing a data validity test based at least in part on the 
power spectra of each range gate bin to determine a 
reference signal for each of the data segments; 
determining a frequency change from the reference signal 
for each of the data segments; 
determining a zero Doppler alignment for each of the data 
segments based on the frequency change; 
determining a power frequency for each of the data seg- 
ments based on the data segment’s zero Doppler align- 
ment; 

determining a Doppler shift f for each line of sight based on 
the line of sight’s respective maximum power fre- 
quency; 

determining a power spectra p for each line of sight based 55 
on the line of sight’s respective maximum power mag- 
nitude; 

determining a line of sight range for each line of sight; 
plotting the line of sight range for each line of sight versus 
the line of sight’s Doppler shift f and power spectra p; 60 
determining east, north, and up wind components for each 
line of sight based on the plot of the line of sight range for 
each line of sight versus the line of sight’s Doppler shift 
f and power spectra p; 

computing a wind speed and a wind direction using the 65 
east, north, and up wind components for each line of 
sight; 


1 

^ n=0 

m = l:M, 


M is a number of range gate bins, n is the data segment 
number, and is the zero Doppler alignment. 

15. The non-transitory processor readable medium of 
claim 14, wherein the stored processor-executable instruc- 
tions are configured to cause the processor to perform opera- 
tions such that: 

determining a Doppler shift f for each line of sight based on 
the line of sight’s respective maximum power frequency 
comprises determining the Doppler shift f as a value 
equal to {fj • ■ ■ wherein 


fm = fo,K + max[K„(/)] - 100 MHz, m = 1 : M, 


^OJC is the reference signal, and F is a nominal value in 
Hz; and 

determining a power spectra p for each line of sight based 
on the line of sight’s respective maximum power fre- 
quency comprises determining the power spectra p as a 
value equal to {p^ . . . p^} wherein 

/)^=max[r^(/)],^F, and m=l \M. 

16. The non-transitory processor readable medium of 
claim 15, wherein the stored processor-executable instruc- 
tions are configured to cause the processor to perform opera- 
tions such that: 

the associated INS/GPS data includes an aircraft speed east 
value (V^^), an aircraft speed north value (V^^y), an 
aircraft speed up value an aircraft roll value, an 

aircraft pitch value, and an aircraft heading value for 
each segment for the two or more lines of sight, and 

determining east, north, and up wind components for each 
line of sight based on the plot of the line of sight range for 
each line of sight versus the line of sight’s Doppler shift 
f and power spectra p comprises determining speed east 
speed north (V^^), and speed up (V jyjj) wind 
components for each hue of sight (d) as vectors defined 
by: 
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tions such that the roll, pitch, yaw, elevation, and azimuth 
offset values are estimated at least in part on a Doppler shift in 
a reference ground return range bin and a translated aircraft 
vector for the reference ground return range bin. 

18. The non-transitory processor readable medium of 
claim 17, wherein the stored processor-executable instruc- 
tions are configured to cause the processor to perform opera- 
tions such that the Doppler Wind LIDAR data includes data 
segments for only two lines of sight which are perpendicular 
to each other, and wherein computing a wind speed and a 
wind direction using the east, north, and up wind components 
for each line of sight comprises: 

adding the wind components for each line of sight such 
that: 


to an angle from N; 

determining an aircraft yaw value as the aircraft head- 
ing value minus the angle from N; 
translating the aircraft vector for each range gate bin 
in reference to an adjusted aircraft axis using the 20 
aircraft speed east value (V^^), the aircraft speed 
north value (V^^), the aircraft speed up value 
(Sau\ aircraft roll value, the aircraft pitch 
value, the aircraft heading value, the aircraft yaw 
value, and roll, pitch, yaw, elevation, and azimuth 25 
offset values; 

estimating wind (W) components in reference to the 
adjusted aircraft axis based on the Doppler shift f 
for each range gate bin and the translated aircraft 
vector for each range gate bin; and 20 

translating the wind (W) components in reference to 
the adjusted aircraft axis into Earth reference wind 
components speed east (V^y^), speed north (V^n), 
and speed up (Vwu). 

17. The non-transitory processor readable medium of 25 
claim 16, wherein the stored processor-executable instruc- 
tions are configured to cause the processor to perform opera- 






